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POST-FLIGHT TH THE "J-7" SYSTEM 
1. INTRODUCTION 

The previous post-flight thermal data analyses of the 
"J" systems were based upon the assumption that the data 
obtained from those temperature sensors attached to the vehicle 
skin were approximately correct. These data contradicted the 
temperature predictions based upon the assumed orbital thermal 
environment and the exterior surface thermal paint mosaic. 
This discrepancy between the analytical and the empirical 
data was resolved in favor of the empirical skin-temperature 
data since (1) the thermal environment can have temporal 
variations whose magnitudes and frequency are unknown and 
(2) the exterior surface paint mosaic is subject to surface 
finish degrading effects of unknown magnitude during the 
launching phase of operation. It is known that there are 
errors in the skin-sensor data due to the dissipation of 
electrical energy within the sensors; this was a major source 
of error for the sensors located within the vehicle. Assuming 
perfect thermal contact between the sensor and the skin, the 
relative magnitude of this error depends upon the relative 
amount of skin sensor electrical heating of the immediate 
adjacent skin area as compared to the total heat exchange 
for that particular area. Since the total heat exchange for 
most skin sensor positions was very large, involving direct 
solar radiation, -.albedo radiation, and earthshine radiation, 
the relative error was assumed to be small. The skin-sensor 
data for those vehicle areas not exposed to large heat-exchange 
rates, that is, those always facing dark. space can be expected 
to be somewhat high* since the electrical heat input for those 
sensors become a significant proportion of the total heat 
‘exchange for those skin positions. 
. The vehicle-skin-temperature~sensor data, tape recorded 
throughout an orbit, were used in, previous flights as the 
thermal boundary conditions to check the thermodynamic con- 


6istency of the internal sensor data. The interior components’ 


iThe Associate Contractor in a verbal communication October 1l, 


1963, stated that during a HATS test of an entire system that 
the skin sensors did not follow the skin temperature as indi- 


cated by calibration the les in the low temperature 
ranges. 2 i il Sais ee 
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temperatures were computed for the thermal environment pre- 
sented by the vehicle skin as described in Reference 1. With 
the adoption by the Associate Contractor of the practice of 
calibrating the temperature sensors for electrical self heating 
as a function of time prior to flight, the thermodynamic con- 
sistency of the sensor data improved considerably. The average 
variation of the internal temperature data ‘from that computed 
for the "J-5" system flight was approximately +1° F (Ref. 1, 
Table II data) while that for the earlier "J-1" flight was 
approximately +9° F (Ref. 1, Table III data). The flight 

data, tape recorded through an orbit, were necessary for this 
type of analysis since the skin temperatures, which vary con- 
siderably during an orbital revolution, are used to define 

an orbital thermal environment for the interior. 

The Vidya Division of Itek has not received any tape 
recorded data for flights "J-9" and "J-7". Following a request 
of August 27, 1964, to the Associate Contractor by Vidya for 
the reduction of the "“J-7" tape recorded flight data, the 
Associate Contractor submitted to the customer a cost estimate 
for that reduction. Vidya was informed September 28, 1964, 
that the request of the Associate Contractor for funding for 
the data reduction was denied. The data received, therefore, 
have been only those acquired at a single orbital position and 
do not permit a thermal analysis of the previous type since 
the entire boundary conditions cannot be established. 

For the "J-7" post-flight thermal analysis, therefore, 
it was necessary to assume an average thermal environment due 
to the solar, albedo, and earthshine radiation heating to 
which the satellite was subjected during operation. The 
temperatures expected for various sensors were computed for 
the assumed thermal environment and for the exterior paint as 
it was designed. This type of computation can be expected 
“to result in larger differences between the measured and the 
computed internal temperatures than those for the computations 


based upon the use of the measured skin temperatures as the 








boundary conditions. This is due to the inclusion of additional 
variables, the unknown variation in the space thermal environ- 
ment and the unknown variation, if any, in the exterior thermal 
paint mosaic. Thus, an assessment of the thermodynamic con- 
sistency of the internal sensor instrumentation data and the 
Getermination of the sources of thermal errors cannot be accom- 
plished as was done for previous flights (Ref. 13 for example). 
These difficulties with the current data will be explained in 
detail in subsequent portions of this report. 

This report is the last of a series of six thermal analysis 
reports foilowing the initial "J" system thermal design report. 
It is appropriate, then, that a résumé be included in this 
report describing the contributions made by the thermal analyses 


performed. This resumé is presented in Section 5. 


2. METHOD OF ANALYSIS 

The data obtained from the Associate Contractor includes 
the "J-7" system post-flight thermal data as reduced from the 
telemetered raw data obtained for a single orbital. position - 
corresponding to Vandenberg Air Force Base, the paint mosaic 
design for the "J-7" system, and the orbital B angle as a 
function of the number of orbital revolutions of the system. 
This material is presented in Appendix A. 

The 130-node mathematical model of the system, together 
with the resulting solar absorptivity and thermal emissivity 
of the skin for the designed paint mosaic, were used to compute 
transiently the system nodal temperatures over an orbital 
period. The orbital 8 angle used for these computations was 
40° corresponding to the system orbital revolution number 46. 
(Appendix A). A detailed description of the computational 
procedure is given in Reference 3. The heat fluxes for the 
thermal environment at an orbital 8 angle of 40° were computed 
by the methods described in Reference 4; for their computation 
an average terrestrial reflectivity (or albedo constant) value 
of 0.39 was used for the near polar orbit of the "J-7" system; 
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the orbital inclination with the equator for this flight 
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was 85° (Appendix A). 

The orbital position reference point for the computations 
was the system's emerging transit of the terminator. In order 
to relate the transiently computed temperatures to the data 
acquired at Vandenberg Air Force Base, it was necessary to 
derive a mathematical expression in terms of the orbital para- 
meters furnished to Vidya for the orbital time required for 
the satellite to traverse the path from the terminator to the 
data acquisition station. This derivation is presented in 


Appendix B. 


3. RESULTS OF THE ANALYSIS 

A comparison of the in-flight measured temperatures for 
the vehicle interior as obtained from the data acquired at 
Vandenberg Air Force Base with the computed nodal temperatures 
for the nodes, which best approximate the system's temperature 
sensor locations, is presented in Table I. The flight data 
utilized for comparison is that for orbits 40 and 47. These 
orbits have 6 angles of 40.5° and 39.9°, respectively. The 
temperature differences resulting from their departure from 
a £ angle of 40° should be less than 1° F in the predictions.? 
The positions on the instruments corresponding to the nodal 
numbers are presented in the Nodal Breakdown Chart of Refer- 
ence l. 

A graphical comparison is shown in Figures 1 and 2 for 
the measured temperatures on the vehicle skin, barrel and 
fairing, and the model computed temperatures for the orbital 


position corresponding to the data acquisition station. 


4. DISCUSSION OF THE RESULTS 
4.1 "J-7" Interior Temperatures 

The data comparison shown in Table I indicates an overall 
average computed interior temperature that is appr&ximately 
a 


“In Reference 1 it is shown that a shift of 53° in 6 angle 


produces about a 30° F shift in mean orbital skin temperature. 
For a linear interpola © F per degree 


of 6 angle shift. 
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5° F lower than the average of the system interior measured 
temperatures. By comparison, the temperature level computed 
for the “J-5" system was within 1° F of the mean of the 
measured temperatures (Ref. 8). This small difference is 
expected for the use of measured skin-temperature values as 
boundary conditions if the interior data is thermodynamically 
consistent with the skin-temperature data and the thermal 
behavior of the shielding and interior surface finishes is 
understood. 

Although this agreement to within 5° F between predicted 
and measured temperatures is quite close considering allowable 
tolerances of +10° F about a prescribed mean value, it should 
not be inferred from these data alone that the mean orbital 
environmental parameters employed or the skin finish character- 
istics utilized will always be appropriate. If this were 
true the agreement to within about 5° F between predictions 
and actual data should be obtained for each flight. For the 
purpose of comparing the above results with that for another 
flight, the flight data for the "“J-9" system which had the 
same external paint mosaic® as "J-7" and a similar orbit* is 
compared in Table II. The data presented there are for orbital 
passes numbers 9 anda 16 having 8 angles of 40.6° and 39.6°, 
respectively. The equations of Appendix B were used to cor- 
relate the orbital position with respect to the terminator 
for the data with the corresponding computational time in orbit. 
The predicted temperatures for the data comparison were within 
1° F of those computed for the "J-7" flight. The computed 
temperatures, however, are on the average about 10° F lower 
than the measured ones. Thus, the "J-9" system average temper- 
ature is about 5° F higher than that for "J-7". The differences 
for the two systems may be due to differing space thermal 
environments caused by earth-cloud cover variations or they 





3Ccommunication from the Associate Contractor. 


*tThe orbital inclinations were 85° and 80° and the perigee 
altitudes were 99 N.M. and 84 N.M. for "J-7" and "J-9", 
respectively. The apogee altitude for "J-7" was 260 N.M.; 
for "J-9" it was 261 N.M. This information was obtained 
i the flight report files of the Itek Field Service at 
A/P. 
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may be due to variations in vehicle system parameters such 

as that caused by degradation of surface finishes, variation 
in instrumentation performance or to data read-out and reduc- 
tion errors. Also, a combination of the above factors may be 
responsible. 

From the data currently obtained from the "J" systems 
and the present state of the knowledge of the space thermal 
environment, the effects of the above two sources, environment 
and hardware, on temperature variations cannot be separated 
even approximately. The possible variations of the system 
mean temperature level due to the temporal variations of 
earth-cloud cover from its mean can be only roughly estimated 
until experimental knowledge is obtained concerning the fre- 
quency and magnitude of these variations. This estimate is 
made in Appendix C and indicates that the “expected” variation 
of system temperature levels for different flights due to 
"expected" variations of cloud cover is about +6° F and -4° F. 

These estimated environmentally induced temperature 
variations might well explain the 5° F difference between 
the mean-temperature levels of “J-7" and "J-9" assuming, of 
course, that the variations actually occurred and are not | 
instrumentation errors. However, there exists an approximate 
10° F difference in level between the theoretical results and 
the measured results” for "J-9". If there were a strong de- 
pendency of the theoretical results upon the absolute value 
of the assumed mean terrestrial reflectivity, or albedo con- 
stant, then a strong possibility would exist for an error in 
the theoretical value of the mean terrestrial reflectivity . 
“The theoretical results were obtained for the use of a mean 

solar constant. The actual value of the solar constant dur- 
ing June is smaller than the mean by 3.3 percent. Thus, the 
actual difference between the theoretical and the measured 
results should be even greater than 10° F since the computa- 
tions using 2 lower heat input to the system would result in 
a lower mean temperature. All future computations by Vidya 


Division of this type will include the value of the solar 
constant existing for that flight date. 








used in the computations. However, it is shown in Appendix D 
that for the "J" systems in an orbit having 6 * 40°, the 


mean orbital temperature of the satellite is almost completely 
independent of the magnitude of the mean or long term value 

of the terrestrial refelctivity; only temporal variations in 
the terrestrial reflectivity cause temperature variations. 


Thus, the hardware parameters or perhaps the data itself must 
also be suspected as contributing to the above difference 
between theoretical and measured temperatures. There can be 

no definite conclusions reached by these arguments, however, 
Since the estimated variations in the mean orbital temperature 
due to temporal cloud-cover variations are merely estimates. 
Therefore, as stated above and now rephrased for emphasis, the 
present dearth of knowledge concerning the frequency and magni- 
tudes of variations in the amount of earth-cloud cover prevents 
the separation of the causes of satellite temperature deviation 
from the designed-for value into those causes due to the 
environment and those due to the actual hardware. The separa- 
tion of these causes requiresthe experimental determination 

of the frequency and magnitude of the variations in the extent 
of cloud cover. If these variations are found to be larger, 
then the separation of these causes requires the monitoring 

of the space thermal environment existing for each flight 
either by the means of devices aboard each "J" system or of 

a separate, relatively long-life satellite, especially designed 
for this one task, having approximately the ‘same orbital incli- 
nation as that of the "J" systems. A stronger incentive for 
the improvement of passive thermal control technigues may 
result if it were shown that the temperature variations from 
the expected values now experienced by the "J" systems are not 
due solely to the variations of the space thermal environment 
from flight to flight. On the other hand, if the variations 


in the space-thermal environment are shown to be the source 


of the thermal control difficulties then any improvement in 
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the temperature control will require the use of active thermal 


control methods. 


4.2 "“J-7" Skin Temperatures 

The graphical comparison of the computed skin tempera- 
tures with the measured skin temperatures for the orbital 
position corresponding to the data acquisition station is 
presented in Figures 1 and 2. 

Qualitatively, these results are about what can be 
expected for the present model and the current hardware sensor 
instrumentation. The computed barrel temperatures for the 
sun side are expected to be somewhat lower than the measured 
values for a daylight pass over the data acquisition station 
for the reason explained in Appendix H. It would follow from 
this same reasoning that the computed temperatures for those 
same nodes should be higher than the measured values during 
a night pass; however, the electrical self heating of the skin 
sensors, discussed in the Introduction, will dampen the rate 
of the indicated temperature drop during the night pass°®; 
consequently, the measured results are probably high. 

For the fairing section, the conical SESEION: the computed 
temperatures are expected to differ from the measured ones, 
but without specifically including the first recovery unit in 
the mathematical model the algebraic sign of the difference 
cannot be predicted. The explanation for this is also presented 
in Appendix H. 

For the cold side of the vehicle, the side not exposed 
to direct solar radiation and for the opposite side as well 
during the night pass, the measured temperatures are higher 
than the computed ones as shown in Figures l and 2. The 
probable reason for this is that given above and, also, pre- 
viously in the Introduction with the supporting statement of 
the Associate Contractor (Footnote 1); namely, the electrical 


self heating of the skin-temperature sensors introduces a 





Spuring the daylight pass this self -heating will accelerate 
the indicated temperature rise, but the total heat exchange 
rate due to the direct solar radiation is of such a greater 
magnitude that the effect of the sensor self heating is 


negligible by comparison. 
erence aes. 
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a positive temperature error in the data for the low tempera- 
ture ranges. 

The above discussed differences between the computed, 
instantaneous, transient temperatures and measured temperatures 
at the single, satellite orbital position shown in Figures 1 
and 2 are of a magnitude of from 0° to 30° F. The ampli- 
tude of the skin's transient temperature over an orbital period 
is from Figure 1 about 80° F; this is the difference between 
measured, instantaneous temperatures at a body angle of 60° 
for the day and night passes for the data acquisition station 
position. From the 35 percent error in temperature amplitude 
expected, discussed in Appendix H, p. H-3, an error of about 
30° F could be expected in the computed temperatures. Also, 
in Appendix H it is shown that this error in the computed 
instantaneous skin temperature results in an error of only 
about 1° F for the interior temperatures. Thus, the errors 
in the computation of the transient, skin temperatures result 
in a negligible error in the computed interior temperatures. 

The effects of the variation .of cloud cover as discussed 
above in Section 4.1 will apply to the exterior transient 


temperatures also. 








5. RESUME OF THE CONTRIBUTIONS TO THE "J" PROGRAM RESULTING 
FROM THE THERMAL ANALYSES 


5.1 Thermal Design 

The thermal design analyses performed by Vidya Division 
in liaison with the Associate Contractor resulted in a passive 
thermal control design for the system that would permit system 
operations at any orbital ¢ angle. With this design, an 
interior temperature between 70 and 80° F should be maintain- 
able for any 6B angle (Ref. 1, Table V). ; 

To attain this performance it is necessary to employ a 
radiation shield liner having a thermal emissivity of about 
0.12. It is required that the mean orbital skin temperatures 
of the vehicle at af = 53° be about 110° F. This was pointed 
out explicitly in the original design report on page 6 of 
Reference 1. Quoting from that page, "... the best instrument 
temperature ranges are experienced with a nominal skin temper- 
ature-of 110° F." The underlined portion is underlined in 
the original text. : 

Unfortunately, the thermal radiation shield fabricated 
for the system, as measured by the Asseciate Contractor, had 
an emissivity between 0.2 and 0.4. The variation of thermal 
performance with the numerical value of this emissivity is 
nonlinear; the mean value of the emissivity, 0.3, is com- 
pletely inadequate to provide the required shielding. Com- 
Putations made by Vidya Division using the 130-node mathemati- 
cal model show that this shield is but little more effective 
in providing the desired thermal performance than a black 
shield would be (Ref. 13, Section 5.1). 

For flights more recent than the "J-7" system flight a 
new radiation shield has been incorporated. This shield 
appears visually to have a lower thermal emissivity than the 
value 0.3 described above. However, for the flight thermal 
results for the first two systems containing this shield as 
reported qualitatively to Vidya Division, it appears that the 


designed mean orbital vehicle skin temperature was not in the 
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required 110° F range at a § = 53°, but was much lower instead. 
The resulting interior temperatures as indicated by the flight 


data (again, a qualitative report to Vidya Division) were low, 


in the 45° F range. This kind of result can be expected to 
be obtained in the future unless the original design philosophy 
of maintaining high mean skin temperatures is followed’. 


5.2 Post-Flight Thermal Data Analyses 

An important contribution to the thermal performance of 
the “J” systems has been the improvement in the accuracy of 
the interior temperature sensor data. 

Originally, the Associate Contractor noted that, during 
ground tests, variations as a function of time occurred in the 
temperatures as indicated by the interior sensors. A test was 
then performed by the Associate Contractor using a sensor 
attached to a metal honeycomb specimen and from the tempera- 
ture-time history data resulting from this test, a calibration 
curve was obtained to correct the flight data (Ref. 11). 

Vidya Division also performed a test of a sensor attached to 

a metal honeycomb specimen having different dimensions from 

that used by the Associate Contractor (the same type of honey- 

comb material was used as that used in the fabrication of the 

instrument main plates). A different temperature-time history 

for the sensor was obtained by Vidya as compared to that 

obtained by the Associate Contractor. A mathematical model 

7In a recent verbal communication from the Associate Contractor, 
information was received that the third "J" system having the 
new low emissivity radiation shield had, in flight, an interior 
temperature of about 70° F. Also the information relative 
to the external paint mosaic was that the sun side of the 
vehicle had an a/e of about the same numerical magnitude 
as for the previous systems (* 1.1) and the shadow side of 
the vehicle on the side not exposed to the sun was entirely 
gold covered. An approximate computation by Vidya Division 
for this thermal mosaic indicates that a mean orbital skin 
temperature of about 105° F should have been obtained at B = 53°. 
Thus, the thermal design for the most currently flown system 
(circa January 1965) appears to be in line with the original 
design philosophy. The question remains, however, of why it 
was necessary to use a trial-and-error method by which the 
two previously flown systems in which low emissivity shields 


were employed ran cold, to obtain the thermal design described 
one and one-half years ago in Reference 1? 
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of the sensor as attached to a metal honeycomb structure was 
used to analyze the test results. This analysis revealed 
that the two tests, Vidya Division's and the Associate 
Contractor's, were eguivalent, the difference between them 
being introduced by the specimen dimensions used. The basic 
information obtained from this analysis was that one single 
sensor calibration could not be used for the system because 
the amount of sensor self heating was dependent upon the 
sensor's position on the instruments. The Associate Contractor 
then began individual sensor calibration in situ during the 
vacuum chamber tests prior to flight. 

The improvement in the accuracy of the interior sensor 
data is illustrated by the previously reported (page 2 of 
this report) decrease in the magnitude of the measured temper- 
ature level variation from that predicted from thermodynamic 
principles. 

A continuing flight to flight monitoring of the interior 
sensor data was conducted by comparing them with the results 
of the 130-node mathematical model. The measured skin sensor 
data were used as boundary conditions to ascertain the con- 
sistency of the internally measured results with thermodynamic 
Principles. For each flight for which the measured skin data 
were obtained the results were compared in tabular form, as 
for example, Tables I and II of this report. The purpose of 
this monitoring was to provide an independent quality check 
upon the measured data and thermal performance of the instru- 
ment. This is indispensible for a new system since there are 
nearly always some initial system flaws that must be corrected. 
The above described temporary correction for sensor self 
heating (a new system of sensor instrumentation is necessary 
for a basic correction) is an example of the use of an inde- 
pendent monitoring of the data for the initial break-in period 
of a new system. Also, this has resulted in an effort to 


obtain continuity from the initial thermal design effort 


through this break-in persoc SEPRET 
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The written reports for this monitoring function are 
usually quite late in time compared to the system flight 
schedules. The contents of the reports are not the entire 
fruits of the work, however, the liaison activity between the 
two Associate Contractors to improve the system constitutes 
a large part of the results. 

Once the flaws in a new system have been revealed the 
flight to flight monitoring of the data is no longer neces- 
Sary.- Thus, Vidya does not recommend a continuation of this 
effort on a routine basis. It is recommended, though, that 
when changes are made in the internal thermal control surfaces 
or in the sensor instrumentation that an independent check of 
the flight data by Vidya be made for the flight of the system 
in which the changes were incorporated. 

Other relatively short-term tasks performed during these 
post-flight data analyses by Vidya have been the following: 

a- An analysis of the thermal effects of air bubbles in 
the cement used to attach the sensors to the instru- 
ment (Ref. 1). (The effects were found to be detri- 
mental to sensor accuracy and a recommendation was 
made to Itek that tighter quality control be employed 
in the attachment of the sensors.) 

b. An analysis of the effects of the transient exposure 
of the lens to the exterior environment during 
operation (Ref. 2). (The resulting temperature 
changes were found to have a negligible effect upon 
record quality as measured by the image shift from 
the record.) 

c. An analysis of the effect of temperature level upon 
record quality (Ref. 2). (This analysis required 
the use of detailed temperature histories, the tape 
recorded flight data, for several flights. The 
detailed thermal data reduction was discontinued 
beginning with the "J-9" system and this analysis 


procedure was aborted.) 








6. CONCLUSIONS | 
It is concluded from the results of the analyses of the 
"J-7" and "J-9" systems that . 

a. The differences between the mean of the measured 
temperatures and the mean temperatures computed 
by Vidya Division for the solar constant's value 
existing at a given flight date for the "J" systems 
are due either to temporal variations in the space 
thermal environment or to variations in the values 
of the hardware design parameters. 

b. In order to evaluate numerically the relative effects 
of the two causes (Item a above) of the variations 
in a satellite's thermal performance, it is neces- 
Sary to determine the magnitude and frequency of 
the variations in the space thermal environment 
and possibly to monitor the space-thermal environment 
during each flight. 

c. Should the necessity arise to improve the accuracy 
of the transient temperature computations for the 
vehicle skin it will be necessary to increase the 
number of nodes in the mathematical model to include 
the skin reinforcing structures as separate nodes. 

From the thermal analyses of the first six-system flights 

it is concluded that® 

a. The temperature sensor instrumentation should be 
changed to eliminate the necessity of individual 


sensor self-heating calibrations’. 





Sat the current date (February 1965) the emissivity of the 
internal radiation shield of the vehicle is on the order of 
that specified in the initial design report, therefore, these 
are conclusions not concerning the shield. 


“In Reference 2, Section 3.3, it is pointed out that the cali- 
brations change as a function of the environmental tempera- 


ture level; hence, calibrations at one level used for system 
temperatures at different levels will result in data errors. 
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b. Air bubbles should be entirely eliminated from the 
temperature sensor cement (Ref. 13). 

c. Greater accuracy in the computation of the internal 
temperature distribution of the system reguires a 
more detailed nodal breakdown (more nodes), experi- 
mentally measured values of the thermal contact 
resistance across riveted and bolted structural 
joints, the experimentally measured values of the 
specular reflectivities of the different surface 
finishes, and the inclusion of the specular nature 
of the surface finishes in the computations?°. 

d. A post-flight thermal data analysis is not necessary 
for each future flight of the "J" systems. Periodi- 
cal analyses should be made however to assure the 
maintenance of quality control in data acquisition 
(instrumentation plus data readout and reduction). 
Also, the results of new sensor instrumentation 
and new interior surface finishes should be checked 
by an analysis of the flight results. 

From the flight results to date (circa February 1965) 
and the indications that a period of one and one-half years 
has been required to attain in the hardware the initially 
specified design properties for the "J" system, it is con- 
cluded that there is considerable room for improvement in 
the quality of the liaison between design groups and opera- 


‘tions groups. 


7. RECOMMENDATIONS 
The recommendations based upon the conclusions of 
Section 6 are as follows: 

a. That the magnitude and frequency of the variations 
in the space-thermal environment of the "J" systems 
be measured experimentally. 

1°n11 of these items need not be accomplished simultaneously 


to improve the accuracy; for example, the increase in the 
number of nodes will in itself increase the accuracy as 


explained in Appendix H of this report. 
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b. That the type of temperature sensor instrumentation 
existing in the systems through "J-7" should be 
replaced with instrumentation that eliminates the 
sensor self-heating error. 

c. That post-flight thermal data analyses be made for 
perhaps only one in ten future flights for the pur- 
pose of maintaining an independent check upon the 
guality of the flight data and thermal design. 

ad. That post-flight thermal data analysis be made for 
the flight of a system in which new temperature 
sensor instrumentation or new internal surface 
finishes are employed. This analysis should be 
performed as soon as possible after such a flight 
in order to permit its results to be applied to 
the design of the immediately succeeding flight. 

e. That post-flight thermal data analyses be performed 
for any flight for which abnormal thermal performance 
is suspected from the thermal data. 





Scent 
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TABLE I.- COMPARISON OF “J-7" FLIGHT DATA WITH COMPUTATION RESULTS. 


DAY _ PASS NIGHT PASS 









System S_ 
Sensor 
1-I-3 59 50° F{ 55° F +5° F 49° F| 56° F +7° F 
1-1-4 42 54 52 -2 59 453 ~6 
1-I-5 61 a 59 +2 59 60 +] 
1-I-6 56 69 63 -6 73 64 -9 
1-I-7 52 65 60 -5 66 60 -6 
1-1-8 61 61 59 -2 63 60 ~-3 
1-I-9 60 67 63 -4 71 63 -8 
l-I-10 52 66 60 -6 64 60 ~4 
1l-I-11 16 66 57 -9 64 57 -7 
1-I-12 59 =e 55 +4 52 56 +4 
1-I-13 LS 68 58 ~10 70 57 -13 
2-I-3 70 76 62 -14 76 62 -~14 
2-1-4 47 72 69 -4 74 60 ~14 
2-I-5 72 65 58 -7 66 59 -7 
2-I-6 67 60 25 -5 61 56 -5 
2-I-7 63 65 57 -8 63 57 -6 
2-1-8 72 58 58 0 67 59 ~8 
2-I-9 71 55 53 ~-2 58 54 -4 
2-I-10 63 66 a) ~9 64 57 -7 
2-I-11 37 a7 58 +1 61 56 -5 
2-I-12 70 72 62 -10 72 62 -10 
2-I-13 36 64 55 -9 66 55 -11 
SP.-l 75 54 59 +5 57 60 +3 
- SP.-2 75 60 59 -1l 63 60 -3 
CL--l j; 7 70 68 ~2 72 72 0 
CL.-2 74 72 68 -4 75 72 -3 
T. measured temperature 
Tu computed, model nodal temperature 
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TABLE II.- Comparison of "J-9" Flight Data with Computation Results. 





DAY PASS NIGHT PASS 


Approx. 
System Node 
Sensor | Number 





_ Ll-I-3 59 
1-I-4 42 

1-I-5 61 

1-1-6 56 75 63 -12 79 64 fF eS 
1-I-7 52 76 60 -16 79 60 -19 
1-I-8 61 69 59 -10 72 60 -12 
1-I-9 60 71 63 -8 76 63 -13 
1-I-10 52 73 60 -13 72 60 -12 
1l-I-11 16 75 57 -18 72 57 -15 
1-I-12 59 57 55 ~2 59 56 -3 
1-I-13 £5 a7 58 -19 79 57 -22 

- 2-JI-3 70 74 62 -12 74 62 ~12 
2-1-4 47 71 69 -2 73 60 -13 
2-I-5 12. 65 58 -7 68 59 -9 
2-I-6 67 61 25 -6 62 56 -6 
2-I-7 63 68 57 -1l 68 57 -11 
2-I-8 72 63 | 58 ~5 64 59 -7 
2-I-9 71 57 53 -4 58 54 -4 
'2-I-10 63 74 57 -17 70 57 -13 
2-I-11 37 58 58 0 60 56 -4 
2-I-12 70 66 62 -4 68 62 -6 
2-I-13}° 36 72 55 -17 73 55 -18 
SP.-1 75 53 59 +6 55 60 -15 
SP.-2 v2 61 59 -2 62 60 -2 
CLe-l 74 75 68 -7 | 71 72 +1 
CL.-2 74 73 68 -5 69 72 +3 

T. measured temperature 
Ty Computed, model nodal temperature 
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APPENDIX A 


THE "J-7" SYSTEM POST-FLIGHT DATA RECEIVED 
FROM THE ASSOCIATE CONTRACTOR 
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APPENDIX B 


DERIVATION OF THE SATELLITE'S ANGULAR POSITION FROM 
THE TERMINATOR PLANE. 


In Figure 1. the origin O is assumed to be at the earth's 


center. In addition set ‘4 


E the north ecliptic pole ° 

N a unit vector along the north equatorial pole 

‘e) a unit vector along the orbital normal 

S a unit vector along the earth-sun line 

P | a4 unit vector along the earth-satellite line 

T a unit vector in the direction of the ascending inter- 


section of the orbit and terminator plane 


¥ 


Vv a unit vector in the direction of the vernal equinox 
Pp’ the projection of P onto the equatorial plane 
Ss’ the projection of S onto the equatorial plane 


‘ 


The following angles are assuncd to be known quantities: 

sun's right ascension 

i orbital inclination with the equator 
the angle from the north pole to the satellite position 
(90 minus its latitude) 

Z 23-1/2° 
(180 - ::) where :.. is the angle between the ethical 
plane normal as defined by the right-hand rule and the 


earth-sun line. .° can vary from 9 to 180° 








A relationship between the known quantities and the angle $6 


is required. The angle 
\ 


is the position of the satellite in its 


Orbit as measured from its ascending cross of the terminator plane 


ae se 0-0. O-P: 
Or Ss Ox P = A A A on 
7 S:+0O $-P 
sin (180 - +') cos (- - 180 + 90) = cos y - 0 
Sin &' sin © = cos ¥ (1) 
i 4 aS nay 
yo nA», ;N-N N: S 
NP Nv S- = i, : 3 
«| Pen P-s. 
sin & sin x cos (" - ~) = cos + - cos &€ cos x { 2) 
4 A n~ ~ | 
a ae N-.N N-S | 
‘NSE Nes Sone: FS a eI 
Y 'E.N E-S ; 
sin Z sin x cos (90 + 7) = 0 - cos Z cos x 
ctn Z 
tan 2 sence (3) 
(. can be either lst or 2nd quadrant) 
4 wae ey ose 'N-N N:-S 
(Nx 0 | (Nx Sto = | . a 
; -O-N Oo -S 


sin i sin . cos (.-90 - °) 


cos (90- + ') 


cos (180 = ?') 


~cos (' - cos i cos x 


Sin i sin x 


- cos icos « 


Hi 


(4) 





B=3 


From right spherical trigonometric relationships the 


following equations are obtained. 


cos Z 


tan = arenas (5) 
csc 360 - ( - -) = e@6¢(a = fy cs <a (6) 
Ss anes ; ctn = 
From these $1x relationships the value of «& can be com- 


puted for the given orbital data. However. there will arise 
ambiguous values for ‘: instead of two values corresponding to 
the sun-lit side and the night side, there will be in general 


four solutions. It is necessary, therefore, to derive equations 


to test the value of for their applicability. These are 
derived as follows. e 
bog oe aN 
; : 0-0 O-P : 
Oo ~ N ; O.-% P = i Mii Pe ah: 
PNee N -P i 


(sin i) (1) cos ( + S + 90 - 180) = cos = = 0 


cos = = sin iocos (7 + <<. = 99) 

cos = : Sin icos {90 -- - ‘) 

cos = = sini sin (+: + A) (7) 
Given :', 1, and -. compute 4 from the law of cosines as 


given in any spherical trigometry text. 








? 


cos \ = -cos :'' cos i+ sin fi’ sin i cos J 


é cos » + cos {' cos i 
ee sin * sin l (8) 
There is only one acceptable value of A. that is, & 7 1807 


Since it is the angle for crossing the terminator as measured from 
the ascending equatorial node. 
Ordinarily the value of 3' is not given as such; instead a 


value <« = .90 - °': dis given. This value for an afternoon 
launched, north-south orbit is such that the value of ‘¢, the 
angle between the orbit plane normal and the earth-sun line, is 
greater than 90° and 2 = + 90° ana accordingly 7' = 180 - 2° 
- 180 -~ 90 - © or =' = 90 - ©. The orbit-plane normal in this 
case points away from the sun direction. For a morning launch, 
again north-to-south. the orbit-plane normal (right-hand rule) 
points toward the sun direction, the angle G «< 90° and hence 

90°. Also. for a morning launch, again north-to-south, the 
angle (.- :- 90') will always be greater than 180°. For an 
afternoon. north-to-south launch. it will be less than 180°. 

The two values obtained for (. - ¢“) for the sin (1 - 0@) 
are both pertinent solutions. One value corresponds to the 
orbital point latitude of (90 - ¢) for the sun-lite side and the 
other to the dark side of the earth. For the angular inclination, 
i ° 90 the smaller value of (. - *) corresponds to the sun-lit 
Side while for i 90 the smaller values correspond to the 


dark side. 





ud 











APPENDIX C 
AN ESTIMATE OF THE “J” SYSTEM TEMPERATURE VARIATION TO BE 
EXPECTED FOR VARIATIONS OF EARTH-CLOUD COVER 

An estimate of the effect of the expected cloud-cover 
variation upon the "J" system mean orbital temperature can be 
obtained by using the results of the measurements of the earth 
albedo by Danjon and Dubois (Ref. 5). Their measurements, 
made over a period of several years, indicate that the earth's 
albedo varies over a range of about +0.12 and -0.08 from its 
mean value. Their measurements were averaged in blocks of 
1 month each and, hence, possibly laroer variations occurred 
for shorter periods of time but in the averaging process they 
are not apparent. The assumption is made here that the magni- 
tudes of the variations in their measurements from the mean 
are indicative of the magnitudes of terrestrial reflectivity 
variations that can be expected to exist from one "J" system 
flight to another. 

The equation derived in Appendix E relating the system 
orbital mean-skin temperature to the environmental heat fluxes 
is : 


x 


pt x ra r= + ~ 
CT pS MG a, 


The albedo flux, G,, incident on the satellite is related to 


A? 


the earth reflected radiation by the equation 


where ‘ 
Lan is a constant derived frem the satellite geometry 
and the particular satellite orbit involved (this 


is apparent from the equation for the albedo flux 


cf Appvendix E.) 





C-2 


a is the terrestrial reflectivity or albedo constant 


E is the solar constant 
For a given satellite orbit, including its * angle’, the 


variation of the value of Gy is just proportional to the 


variation of a. The values of G. and Gy as a function of 


® angle for an earth oriented cylindrical satellite af an 


altitude of 150 statute miles are presented in Figures 1 and 


2. For a & angle cf 40° and for a = 0.39 


r- 


Ce 16.4 Btu/hr- ft? 


G. 69.3 Btu/hr-ft* 


The use of the earthshine equation of Appendix F and the use 
of the equivalent earth black body temperature of 452° R 
results in 

G, = 28.4 Btu/hr-ft® 
The "J-7" system designed paint mosaic for the cylindrical 


barrel had an w/c = 1.086. The mean orbital skin temperature 


computed by the above equation results in 


(= % 
r | = 516° R or 56° P 
Lad 


The increase of a by +0.12 in accordance with the 


positive variation found by Danjon and Dubois results in a 


e ' = O.51 = — 2 
value of CG A of G A (6-39) (16.4) 21.6 Btu/hr- ft for 


which the use of the above equation results in 


= | = 522° R or 62°F 


This is an increase cf 6° F in the mean temperature. 





‘The orbital t angle is the angle subtended by the orbital 
plane with the earth-sun line. 





C-3 


For the negative variation of -0.08 found by Danjon and 
Dubois in their measurements of the albedo constant similar 


computations to those above result in a value of 


=> by ; 
T° = 512° R or 52° F 


“This is a decrease by 4° F. 

For the above computations the assumption is made that 
the earthshine flux remains constant during the terrestrial 
reflectivity, or albedo constant, fluxuations. Since the 
earthshine is the solar energy initially absorbed by the 
earth and then reradiated back to space the energy balance 
requires that the absorbed solar energy is proportional to 
(1 - a) where a is the wean terrestrial reflectivity, or 
albedo constant. For a steady-state condition this is true; 
however, for short time period variations of a (a few 
days) the thermal inertial mass of the earth is such that the 
earthshine flux is relatively constant. Albrecht Computes a 
heat reservoir of about 1.6x10° Btu/hr per square foot of 
the earth's surface (Ref. 6, p. 2-2). Lettau and Haugen aig 
the same reference, p. 2-l, state that 100 days of radiation 
from the earth would be required to deplete this reservoir if 
all Seeray inputs to the earth were to cease. Thus, the 
assumption of a constant earthshine flux for a period of a 
few days during a change in the albedo constant of +0.12 or 
-0.08 appears to be valid. 

For the assumption that the magnitude of the monthly 


variation in the earth albedo constant as measured by Danjon 





c-4 


and Dubois represent the variations that might be expected 
to occur over short intervals of time, the resulting mean 
Orbital temperature variation to be expected for the "J" 
system satellite from this source alone, is approximately its 


mean temperature +6°, -4° FP. 
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APPENDIX D 
THE EFFECT OF THE MEAN TERRESTRIAL REFLECTIVITY UPON 


THE ORBITAL MEAN TEMPERATURE LEVEL OF AN 
EARTH ORBITING SATELLITE 


The earthshine flux, Gos to a satellite is directly propor- 
tional to the earth’s emissive power, oT", (Appendix E) where 
T is the black-body temperature of the earth and o is the 


Stephan-Boltzman radiation constant. This can be expressed 


mathematically as 
mn =p 4 
te.” "6 (oT, ) 
where kK. is a constant of the geometry and particular satellite 
orbit. For an equilibrium or steady-state condition the solar 
energy absorbed by the earth is reradiated to space; this is 
expressed by 
. N P 
(47R*) (oT, *} = (TR°)(E)(1 - a) 
where 
R earth’s radius 


E solar constant 


a the mean terrestrial reflectivity (for this case, the 
entire earth, it’is the planetary albedo) 


The left-hand side of the above equation is the total energy 
radiated to space; the right-hand side is the total solar energy 


absorbed by the earth. Accordingly, the earth’s emissive power, 


_ E({l - a 








Thus, 


e 


e 


The albedo flux, G to a satellite is directly proportional 


A’ 
to the solar energy reflected from the earth's surface (Appendix E). 


Mathematically, this is expressed as 


Ga = k ,aE 


where k, is a constant similar to ko The value of ae for the 
albedo flux is more dependent upon the local portion of the earth's 
surface viewed by the satellite than is the a for the earth's 
emissive power above. However, their values for polar orbits are 
close numerically and will be assumed to be equal for this analysis. 

From Appendix E the equation for the orbital mean temperature 
of a satellite is 

oT* = G+ as G, + “Ss G 


e A S 


It is assumed that the emissivity, «, is for the range of values 
encountered for T equal to the absorptivity value. By assuming 


Enae a, = ao, 


ee a a 
a cs, +(2)a, +(2)c, 


Upon substitution of the above expressions for G. and Ga the 


equation becomes 
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Allowing the constants Ky and k, to absorb the constants E/4 


and (a/e)E the equation further reduces to 


¢ 


4 ' f ‘ Qa” 
JT” = kha - a) + ka +(2)G. 


Partially differentiating with respect to 4 results in 


er ee | 
ee oa) ear a 


va e A 


which for Ke = kn vanishes. 


The value of k., defined as k, E/4, is 





For a value of a of 0.39 the value of G., for the earthshine 


equation of Appendix F is G. = 26.2 Btu/hr-ft~ and 
ki = S72 = 43.0 Btu/hr-ft™ 


where ke is a constant of the satellite and orbit geometry and 
is independent of the value of a used. The value of G. varies 
directly as (1 - a). 


The value of ky with the absorbed (a/c)(E) is 


mae (2) &) 


which, for ; = 40° and a = 0.39, G, = 16.4 Btu/hr-ft~ (Appendix C, 
Fig. C.2). And, for ua/e = 1.086, 
Pe 16.4\_ eee 
kn = (1.086) ( 39 )7 45.7 Btu/hr-ft 


SECRET 








and, thus, 


A 


= + 2.7 Btu/hr-ft* 


Qu 


a 


The value of oT at (3 = 40° is 26.2 Btu/hr-ft* for the 


above assumed value of a= 0.39; and for Ge = 69.3 Btu/hr-ft* 


from Appendix C, Figure C.1, and a@/e = 1.086 for the present 


example. 


ot* 


26.2 + (1.086) (16.4 + 69.3) 


119.2 Btu/hr-ft* 


Fake = 514° R 


Replacing the differentials by finite increment symbols,* 


o(sT ) 2 oleT) 2 420.7) Btu/hr-£t? 
da oa 


Values for a given in the literature range from 0.34 to 
0.43 (Ref. 8, p. 166). Hence, an assumption that Aa * +0.1 is 


quite conservative, and, accordingly 


A(oT*) = +0.27 Btu/hr-ft* 


For the above value of <cT* = 119.2 Btu/hr-ft<, the value of 


f— 1/= 5 
a Tt’ os, 1 F 





‘This is correct mathematically for this case since oT* varies 
linearly as a function a. 
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and, hence, the orbital mean temperature of the system is effectively 
independent of the mean or long-term terrestrial reflectivity. 


This is the reflectivity as averaged over a period of several 


weeks. 








APPENDIX E 


SATELLITE, ORBITAL MEAN TEMPERATURE COMPUTATIONS 


The rate of change of the space averaged skin temperature 


of a satellite space vehicle, such as the "J" system, is 


mc gq =a, G +a, G, +e G, - ect’ (1) 
where 
m unit area skin mass 
c specific heat of the material 
Ge direct solar heat fluxes averaged over the skin 
Ga albedo heat fluxes averaged over the skin 
G., earthshine heat fluxes averaged over the skin 
= time 
a. Solar energy absorptivitive 
oan albedo energy absorptivitie 


iT. 


emissivity 
The emissivity, ¢, is applicable for both the absorptivity 
for the earthshine and the emissivity of the vehicle skin 
Since c« for the paint materials and the gold surface finish 
of the "J" system as a function of temperature changes very 
little over the range of the earth's effective temperature, 
454° R, sna-ene maximum skin temperature, about 700° R. This 
is evident from the curves of reflectivity for these materials 
presented in Reference 7. : 

The rate of change of G. and G with the orbit 


A 


precession, about 24° per day, is such that the mean skin 


o APP 





E-2 


temperature averaged over one orbital revolution changes but 
slightly from Bac orbital revolution to the succeeding one. 
There are about 16 orbital revolutions per day such that the 
orbital £8 angles for two successive revolutions differ by 


less than 0.2°. From the curves of G. and G as a function 


A 
of ff angle, Appendix C, Figures 1 and 2, this is seen to 
result-in a change of G. and Ga at ;~ angle = 53° of less 
than two-tenths of one percent; this is considered negligible 
since the temperature varies as the fourth root of the heat 
flux; a 0.2 percent variation in the heat flux at a, 530° R 
mean temperature would cause a temperature variation of about 
0.2° R. The "J" systems have a magnesium skin of approxi- 
mately 0.070-inch thick. The time constant for the skin 
havina aif: * 1 is about 6 minutes, Appendix G. Since the 
orbital period is more than ten times? this, the orbital mean 
temperature of the vehicle skin will be very nearly the mean 
equilibrium temperature corresponding to that orbit's mean 
thermal environment. Thus, the integration with time of 
Equation (1) over an orbital period will result in a mathemati- 
cal reiationship between the vehicle skin mean temperature 

and the mean thermal environment that is applicable to any 


one single orbital revolution chosen from a succession of 


such revolutions. The resulting equation is, since 
period 


as: = zero 








oc * os: re + + 
oe te G. 2» Ga € Ge (2) 


‘The orbital period for the "J" system is about 90 minutes, 
Appendix A. 





m4 
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In general, where Mor Ans and € vary over the venuere sur- 
face, it is necessary to perform the time and space averaging 
as indicated by the bars in Equation (2); in the present case, 
however, the vehicle skin paint pattern is designed to provide 


a constant a_, and ¢€ asa function of position around 


q 
Ss A’ 


the body; this is true to within the dimensions of the paint 


pattern strippingwidth of a few inches. Thus, the a and e¢ 


values can be separated from the terms T*, om and me and 
for the assumption that De = Aa» this results in the equation 


eT = G/e) (6. 1 Ge G. (3) 


The values of 1 and € depend upon the relative area pro- 4 
portions covered by the various surface finishes making up the 
thermal control surface mosaic of the skin. For the present 
systems two surface finishes are used, white aad and gold. 
Let x be the proportion of the skin area covered by the 


Paint having an absorptivity, to? and an emissivity, Ep? and 


(l-x) be that proportion covered by gold having a5 and Eg’ 


Then 
em x tr + (1 - x) a ‘ 
= © =< + (1 - 
- x p ( x) Eg 


and upon solvina for x, 


Meee ee ee ae (5) 





tet 
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For any desired a/e Equation (5) can be used to compute 
the percent of paint (x x 100) necessary for its composition. 
The individual values of a and € can be obtained then from 


the numerator and the denominator respectively of Equation (4). 





APPENDIX F 


THERMAL RADIATION FLUX TO A SATELLITE 


This Appendix consists of a number of Appendices 
originally compiled for Vidya Report No. 51. It is a self- 
contained unit; all references made therein pertain to 


“sections internal to Appendix F. 
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APPENDIX A 


SYMBOLS USED IN APPENDICES B, C, AND D 


area Of orbiting surface 

solar constant, 443 Btu/hr-ft® 
radiation flux, Btu/hr-£t® 
absolute temperature, OR 
albedo of earth, average = 0.35 
area element of earth's surface 
satellite altitude 

absorption coefficient 


mean radius of earth, 3960 statute miles 


orthogonal coordinates 


angle made by the normal to the earth's surface 
element with the line connecting that element 
to the orbiting surface 


angle made by the normal to the insolated earth surface 
element with the line connecting that element to 
the sun 
angle made by the x-y plane projection of the 
normal to the irradiated surface with the 
orbit path of the satellite 
surface emissivity 


angle in the x-y plane made by the plane of the 
orbit to the y axis 


angle made by the normal to the orbiting surface 
with the line connecting it to the earth surface 
element 


angle made by the sun to earth line with the 
z axis : 


Spherical coordinate angles 


0. 
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Angle made by the line of apsides with the earth- 
sun line. For a circular orbit, the angle made 
made by the earth-sun line with the orbit planes 

angle of orbiting surface normal to Zz axis 


the @distance between the earth surface element 
and the orbiting surface 


Stephen-Boltzmann constant 


angle subtended in the x-y plane by the irradiated 
surface normal and the y axis 


angle made by the line of apsides with the Zz axis 


angle made by the normal to the orbit plane with 
the earth-sun line 
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APPENDIX B’ 


HEAT FLUX DUE TO THE SUN'S REFLECTION FROM THE EARTH 


Basic Equation 


The radiant energy from the sun is partially absorbed and 
partially reflected from the earth's surface. The reflected energy 
from any earth surface element, dS, produces a radiant heat flux 
on a satellite which is in view of the reflecting surface element. 


The basic equation for energy reflecting from an earth surface 
element, dS, and striking the satellite surface element, dA, is 


dq = ag £98 4 cos B cos 7 GA ds 


Tp* 


Symbols are defined in Appendix A. 


The satellite coordinates in space with respect to the sun 
and to the earth's surface must be designated. Also, for a given 
satellite-earth-sun geometry the orientation of the satellite sur- 
face point of interest (for heat flux calculation) must be desig- 
nated. 


The position of the satellite will be designated by the 
orientation of the orbit plane in space relative to rhe sun, the 
eccentricity of the orbit, and orbit perigee. 


The surface orientation is designated by the relationship of 
the satellite surface normal angles made in the coordinate system 
to the earth's surface and to the sun in the same coordinate system. 


The position of the satellite determines the coordinate system 
orientation with respect to the satellite-earth-sun system. The | 
earth's center is the origin of the coordinate system; the satellite 
is always on the z axis. The y-z plane is established as the 
plane containing the satellite, earth, and sun. 


In Figure B.] it is seen that the angles 2, A, and y 
determine the satellite position relative to the earth and the sun. 
For a satellite time in orbit of approximately 4 days, the varia- 
tion of the orbit position due to precession or perigee tumbling 
is assumed negligible. Thus 2 and A will be orbit constants 
and wy is the positioning variable. 


The variables ¢ and @ establish the position of any earth 
surface element. 
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/ 
The term cos a is derived from Figure B.2 


Sin @ 
p * sina (h + ro) 
and 
xr sina ; 
—2—__—— = sin(a - @) = sina cos @ - cos a sin @ 
h + ro 
sina |cos 6 + i = cos a sin @ 
1 + (h/r,) 
so 
[2 + (h/r )| sin 6 
tana * = 
1 + (h ro cos - 1 
setting 
y = h/r, 
Rea: a 1 + y)cos 6 - 1 


[qa + y)? - 2(1 + y)cos 0 + 1] 


The term cos B is also derived from Figure B.2. In Figure B.2, 
let (at the earthcenter) 


— 
Zz 


= aounit vector along the z axis 
=> 
s = a unit vector pointing toward the sun 


—)> 
i = 2 unit vector pointing to the surface element of 


earth, dS 
—> —- _ Smad 
2.2 2.85 
(z xr.) . (zx 8S) = . 
° — => — <> 
ro: 2 ro: Ss 


sin 9 cos ¢ sin O. = cos B ~ cos @ cos 9. 
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hence 


cos 6 ™ sin 9 sin O5 cos ¢ + cos 6 cos 95 


The term cos jn is also derived from the same figure. 
In Figure B.2, let (at the orbiting surface, aA) 


—> 


z * a unit vector along the z axis 


n ™ a unit vector along the normal to the back of the 
orbiting surface 
c = a unit vector lying along p pointing toward 


element dS 
—_> => <p 
z z ae «| 
—_> = —_> = 
(z xn). (z xc) ® 
—_ =o» —_-> = 
n. Zz m.c 


sin & cos | (1/2 - ¢) = (1/2 + x)| sin y ™ cos(r - n) = cos € cos y 


and | 
yeunre- (a - 8) = (T+ @) -a 
From the law of sines is derived 
~ (1 + y)sin @ 
Ne (1 + y)cos 6 - 1 
so 
= o -1 {1 + y)sin 9 
Y eae nen Ee + y)cos 6 - :| 
and 


(1 +.y)cos 6-1 


ay (1_+ y)sin 6 
cos(7 + 6)sin {ean Se: + y)cos 6-1 | 


sin 9 


V(1 + y)? - 2(1 + y)cos @ +1 


sin y = sin(r + 9)cos { ean” [tesa] } 


siny = 





é 
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ae y= = 1 + y) - cos 8 
Va + y)? - 2(1 + y)cos 6+ 1 


cos(r - n) = - cos n ™ - sin y sin € cos(¢ - x) + cos y cos € 


Now 


so 
cos 7 = E 1+ y) = cos 6) cos §& + sin 8 sin ¢ cos(¢ - 
La + y)? -~ 2(1 + y)eos 6 + 1| * 
For 


as = he sin 96 ag dé 


the equation to be integrated is 








G6 ¢ 
TF. Cc 
— 7 san & = ((sin @ sin 6, cos ¢ 
= [(1 + y)? - 2(1 + y)cos e+ 1] 


+ cos @ cos 65) [(2 + y)cos 6 - 1] {c + y) = cos 6] cos & 


+ sin 8 sin € cos(@¢ - »}) ag dé 


Circular Orbits 


For a noon-launched, circular, polar orbit, y is a constant 
and the angle y is just a function of the satellite surface point 
orientation with respect to the orbit path. The y axis in this 
case is tangent to the orbit path. 


The equation is readily integrated with respect to the variable 
p. 








8 
. 
THe a a i _[(1 + y)cos @ - 1) sin @ (+10 6 sin 65 sin(¢) cos € 


= [(a + y)? - 2(1 + y)cos @ + 1)? 


[ca + y) - cos @] + sin® @ sin 6. sin € {[¢ sin 2(9)] cos -x 
ro sin® ¢ sin x + [ca + y) - cos a | cos @ cos 6, cos & 
2 Ss 


$ (8) 
+ cos 6 cos 6, sin @ sin € sin(¢ - x) ad 


Integration with respect to 6 requires, in general, numer- 
ical methods. 


For circular orbits the angle 6, is a measure of satellite 
position and it is convenient to integrate numerically the equa- 
tion for various values of 4g. 


For other than a noon launching the angle yx is not constant 
but is dependent upon position. It is then necessary to divide 
the angle yx into two parts. Looking downward from a position 
above the satellite the following sketch is applicable: 








satellite 


ORBIT PATH 


surface 
point of 
interest 
normal to surface as projected 
into x-y plane 
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In the equation for the reflected flux to the point P the 
angle yx is found in the expression cos(¢ - x) for yx measured 
in the same direction as $¢. 


Since 
y=O+6 
then 
cos(¢? - x) = cos(@ - f - 4) 


For any one satellite surface point the angle 6 is a con- 
stant (for self-orienting satellites) and ¢€ is a function of the 
orbit constants and the satellite position in the orbit. 


By inspection of Figure B.1 it is seen that 


cos 2 


sin ¢ “Sin 6. oe 


This expression also shows the limiting values of the angle 65, 
for an off noon launching. 
Elliptic Orbits 


For an elliptic orbit having the orbit constants e, A, and 
Q the value of € is found as follows: 


In Figure B.1, let 


/ 


n{ 


= aunit vector lying along axis z 


a} 


= a unit vector lying along the line connecting the 
earth center with the sun center 


—> 
Pp = a unit vector lying along the line of apsides 
pointing toward perigee 
then 
og oma = — 
Z. 2 z.p 
(zx). (2 xp) = 
—_ ~ => bad _—> 
Cc. Zz Cc . p 








sin 85 sin y cos € ™ cos \ - cos 65 cos vy 


cos A = cos O5 cos ¥ 


ee ginee sing 


This expression is set equal to 


[oo eee ey. 

cos € = 4/1 ~sin® ¢ = 1 Sea 
sin 

Ss 


and solved for cos O° 


cos 8, = cos y cos } - sin y sin*® 2 - cos* 


In computing the heat flux it is convenient to obtain 4g 
with the above equation from the required position, y, and the 
orbit constants A and ©. The value of sin 95 and cos 9g 
may then be used in the flux equation. For machine computation 
the above expression used to obtain sin 6g from 


1 - cos® Oo, 


does not distinguish among the quadrants and is therefore diffi- 
cult to program. 


The variation of satellite altitude for an elliptic orbit 
requires that y be a variable. 


y =h/r 


°o 


The equation for the radius vector of an ellipse is 


afl - e® 


t“ T+ @ cos y 


The altitude of the satellite is 


i a. oe eee): 
Rone = EG = ee cos y 7 *o 
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and 
é 2 . 
oe h a(l - e*) - es 
Se Se con a 

ty ro (Ql + e cos ¥ 

and 


: B-«—au—e) 
(1 + y) a ro ri + e cos ¥ 


Substitution of this expression into the flux equation for 
(1 + y) extends the equation to elliptic orbits for the assumed 
conditions. 


Integration Limits of @¢ 


The variable @$ is completely independent of the variable 6 
for only that portion of the earth's surface symmetrical about the 
Zz axis which may be free from the effects of shielding or of sun 
shadow. 


If the tangent plane to a point on the curved surface of a 
satellite passes through the earth then all points on the earth 
surface to the back side of the plane are shielded from the view 
from that point. Radiant energy reflected from the earth from 
surface points to the rear of such a tangent plane will not fall 
on that satellite surface point. The line of demarcation on the 
earth surface is that for which the angle n = 90° (Fig. B.2). 


Thus these limits of ¢ can be obtained by setting the 
expression for cos 7 = 0., This results in such limits being 


oc +5 t (cos | PAS ah gees Aaetnt)) 


Sin @ 


The sun shadow line on the earth's surface is the daylight- 
dark line. The @¢ limits for this line are obtained by setting 
the expression 


cos B = 0 
This results in the limits 


= - =i 
o=+ [ 7 cos (ctn 6 ctn 64) | 








f 
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Integration Limits of @ 


The limits required for the angle 9 are those for the inter- 
section of the sun shadow with the y-z plane great circle on the 
earth's surface, the minimum value of 9 for the line of the tan- 
gent plane to the satellite surface point intersection with the 
earth surface, and the limit of the earth surface visible from the 
satellite (i.e., the horizon). 


The point of intersection of the sun shadow with the y-z 
plane great circle of the earth is 


= 1/2 - 6, 


The point of the minimum value of 6 for the tangent plane 
intersection with the earth surface is obtained by equating the 
limits of @ pertaining to the line of intersection of the tan- 
gent plane with the earth surface. 


rae bx (coo +{ Ls vi cog seme) ,) 


sin @ 
thus 
-1 Wen + y) ~ cos e| ctn € 
C +65 + (cos ( ns - r) 
a a “2 1 + y) - cos @) ctn 
g +6 (cos ( sin @ }- +) 
and . 


2 1 + etn= + 1 - y ctn® 2 + 
1 + ctn® ¢ 


@ = cos 


The limit of © for the earth horizon is 


= —~1 Ree ees 
re] cos (x = 3) 
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Figure B-1.- Satellite-earth-sun coordinate system. 
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Figure B-2.- Earth-orbiting surface geometry. 
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Figure B-3.- Direct insolation geometry. 
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APPENDIX ¢ 


HEAT FLUX TO A SATELLITE SURFACE DUE TO EARTH RADIATION 


If the net heat transfer to the earth is zero, that is, thermal 
equilibrium exists, then the total radiant energy from the sun is 
equal to the reflected radiant energy from the earth plus the energy 
radiating from the earth due to its temperature. For these assump- 
tions the earth radiates as a black body having a temperature of 
approximately 454° R. This corresponds to an average radiant flux 
from the earth's surface of about 73 Btu/hr-ft*. 


Referring to Figure B/1 and using the Stephen-Boltzmann radia- 


tion law, the energy striking an element, dA, from the earth 
Surface element, aS, is found from 


eoT* cos a cos 4 
aq alle m ds dA 


p 
where 
cos a. 1 + yieos 9 -] 
p* Lc + y)* -2(1 + y)eos 6 + 1| . oe 
cos 7n *® i+ -— cos a cos € + sin 6 sin € cos(¢ - € - 5 


(1+ y)? -2(1 + yleos 0 +1}%? 
as = : sin @ dO d¢ 


The flux at a satellite surface point is 


@ ¢ 
4 f 4 . 
i ribet uison ta tlsin so, {fa + x - coe @] cos é 


iG + y)* -2(1 + y)eos 6+1 


+ sin 8 sin —€ cos(@ - ¢ - | d¢@ dé 


Integration with respect to ¢ produces 


8 
oo” Pic-aow eda 6a athe 
is J (1 + y)? - 2(1 + y)cos G41" {([a +» ae ie ale 


." 
+ sin 6 sin € sin(@ - € - 6) d 
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C=2 


For the special case of a horizontal flat plate where 
—é = 0° 
the above equation may be further integrated mathematically with 
respect to the variable 6. In general, however, integration with 
respect to @ requires numerical methods. 
Integration Limits of $@¢ 


The only limits imposed on ¢ are those due to the shielding 
effect of the satellite body itself as discussed in Appendix B. 


The limits, as derived in Appendix B, are 


g=r+4t Gorn{ Kae - cos a] ctn t} 2) 


Sin @ 


Integration Limits of 96 


The limit imposed on the angle 9 is that corresponding to 
the earth horizon as viewed from the satellite. 


age a! ey Gene! SA 
@ cos ( ar =) 


Elliptic Orbit 


For an elliptic orbit the altitude varies with respect to the 
satellite orbit position. The substitution 


a(l - e? 


oe ye ro( + cos y) 


as derived in Appendix B must be made in the flux equation. 








APPENDIX D 


DIRECT INSOLATION OF THE SATELLITE 


The radiant energy flux due to direct insolation on an orbit- 
ing surface is given by 


F = E cos 8' 


The term cos B' is derived from Figure B.3 as follows: 


Z = a unit vector lying along the z axis 


S = a unit vector lying along a line joining the 
surface point to the center of the sun 


N = a unit vector normal to the insolated surface 


_ — — = 
Zz Zz Ze 


—> <> — << 
(zx S) . (zx N) = 
— = > => 
S z S.N 


sin O5 sin € cos(é5 + €) = cos B'- cos 95 cos & 


Therefore, 


cos B'= sin 95 sin € cos(6 + €) + cos 95 cos & 
For an elliptic orbit , 
sin 9 = cos y cos A - sin vy Vsin®Q - cos*) 
as derived in Appendix B. 
= cos 2.) 
C sin (ee) 8 


as derived in Appendix B. 


The angles € and 6 are derived from the particular satellite 
geometry for the surface point of interest. Thus, 


FRE [sin 0, sin € cos(6 + €) + cos 6, cos é| 
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There is always direct insolation on some part of the satellite 
except when the satellite is in the shadow of the earth. The limit- 
ing values of 6g for the earth shadow are 


° -1 1 ) a 1 
90° + cos (AS) < 9, < 270 - cos (x45) 


as shown in the following sketch: 














APPENDIX G 


THE "J" SYSTEM SKIN TIME CONSTANT FOR RADIATION HEATING 


The mean rate, averaged over the skin surface, of heat 


gain for the satellite skin is described by the differential 


equation 
mc “a =a. G. +a, G, + € G, - coT* 
where 
m the mass of a unit area of the satellite skin 
C the specific heat of the skin material 
Oe the solar abosrptivity of the skin 
oy the albedo absorptivity of the skin 
€ the earthshine absorptivity and emissivity of 
the skin’ 
G. the instantaneous solar flux averaged over the 
skin surface 
Gy the instantaneous albedo flux averaged over the 
skin surface 
G, the earthshine flux averaged over the skin surface 
C the Stephan-Boltzman radiation constant, 
0.1713x10 * Btu/hr- £t?-°R*-ster 
T the instantaneous mean skin temperature 


* time 
This equation cannot be integrated mathematically for 
a. G. and a, Gp being time dependent functions; ¢ G, 


is not time dependent. In order to obtain a time constant; 





“This equality of the absorptivity and emissivity is explained 
in Appendix E. 
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that is, the time required for the skin temperature Aenee to 
rise or fall to within (1 - 2) of its Steady -etace value 
the orbital mean values of G. and Ga will be used and the 
time required for the skin temperature to approach to within 
(1 - Ay of its steady-state value will be computed. This 
period of time will be defined as the skin time constant. 
Let 
st = (x. G_ + 4S 


+ 
A E G,)/e Cc 


= a constant 
The bars over Gs and Ga denote the use of the mean orbital 


heat fluxes. Then, the differential equation becomes 


3 
nN 





aT . <4 _ at 
Se ae? T 


Integration with respect to © for s* » T* (that is, a 


heating situation) results in the equation 


oa \ (s + T >) (s - T) S/ Ss / 
where 

Te the initial skin temperature 

T the skin temperature at time @ 
For a £ = 40° orbit the values used for the mean orbital 


heat fluxes averaged over space and time to a horizontal 
cylinder in a circular orbit having an altitude of 150 statute 
miles will be those presented in Appendix C, namely, 


69.3 Btu/hr-ft* 


OQ) 
i] 


16.4 


?) 
a 


G_ = 28.4 
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and for a/€ = 1.086 (as for "J-7") and the assumption that 


gee 1.086 (69.3 + 16.4) + 28.4 = 709x10° °R¢ 


0.1713™10 ® 
This corresponds to an equilibrium bonperature os s = 516° R. 
The "J" systems have a mean-skin thickness of 0.070 inches; 
the skin material is magnesium having a density of 109 lbs/ft?; 


the specific heat is 0.25 Btu/hr-lb-°F at approximately 70” <C 


(Ref. 9, p. 291 and 156). 


Assuming an initial temperature of To 450° R, the 


time required for the temperature to attain 


i 


=) (516° R - 450° R) 492° R 


450° R + (1 - 


is from the above equation 


3S = 0.126 hours or 7.6 minutes 
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APPENDIX H 


A QUALITATIVE: EXPLANATION FOR THE DIFFERENCES FOUND TO EXIST 
BETWEEN THE MEASURED TRANSIENT SKIN TEMPERATURES 
AND THOSE COMPUTED USING A 130-NODE 
MATHEMATICAL MODEL OF THE SYSTEM 
Temporal variations in the amount of earth-cloud cover 

cannot be used as a logical explanation for all the differences 
found to exist between the transiently computed temperatures 
of the “J-7" system and those measured. The ditcerenees found 
_in the thermal analysis will be assumed, here, to be actual 
differences despite probable errors in the measured temperatures 
due to é€rrors in the terperature sensor instrumentation and 
data read-out. Temporal variation in the extent of cloud cover 
should result in approximately the same difference between 
the computed and measured temperatures for both the barrel 
section and the conical fairing section of the system. However, 
for the hot side, the sun side, of the vehicle the actual or 


measured differences are of opposite sign, algebraically. 


The Barrel Section 

The accurate computation of the transient temperature of 
a portion of the vehicle skin is dependent upon the correct 
modeling of that portion's thermal mass (the product of the 
physical mass and the specific heat of the material) as weil 
as the correct representation of the surface finish and the 
thermal environment. The mathematical relationship is shown 


in Appendix E. 








The masses of the skin nodes for the 130-node mathemati- 
cal model of the "J" system were obtained by dividing the 
total mass of the vehicle skin plus the mass of the structural 
reinforcing component proportionally according to the skin 
nodal areas. This is approximately correct if the thermal 
conductance between the skin and these structural elements 
is sufficiently large to make the conductive heat-flow rate 
per unit temperature difference between them large compared 
to the radiation heat exchange between the space environment 
and the skin surface. For such a condition the entire nodal 
mass has approximately the same temperature throughout and 
the component elements of a node react as a unit. The heat 
input to the skin surface is effectively distributed through 
the entire nodal mass. This is the assumption made for the 
skin node modeling of the "J" system. Actually, however, 
there exists surface contact thermal resistance across the 
riveted and bolted joints between the structural reinforcing 
elements and the skin. As a result the actual combination of 
skin and structural elements do not react as a single thermal 
mass; the heat content variation of the node is largely con- 
fined to that portion of the node containing the exterior 
skin surface. A more accurate model would include the re- 
inforcing, structures as separate nodes with the appropriate 
thermal conductance between them and the adjacent skin nodes. 


In the equation for the rate of temperature change 


(Appendix E) the rate is inversely proportional to the thermal 
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mass of the re The larger the mass the smaller is the 
rate of nodal temperature change. When the assumed thermal 
mass for computing purposes is larger than that actually 
existing in the hardware the computed rate of nodal tempera- 
ture change will be less than that which actually exists. 
Therefore, for a given radiation environment which varies 
periodically with time the computed nodal temperature will 
always lag the measured temperatures. During the heating 
portion of the cycle the measured temperatures will be higher 
than those computed; during the cooling portion of the cycle 
the measured temperatures will be lower than those computed. 
An excellent example of this is presented in Figure 6 of 
Reference 2. 

The effects of these errors in the transient temperatures 
of the barrel skin upon interior temperatures are quite 
small. The mean skin terperature which controls the mean 
interior temperature does not depend upon the numerjcal values 
of the skin-node masses as is; this is shown by the equations of 
Appendix E. Also, the amplitudes of the temperatures of the 
instrument main plates are approximately 3° F (Ref. 2, Fig. 
9); a 100 percent error in the amplitude of the paves skin 
temperature could cause no more than a 3° F error in the 
anpiatvude of the temperatures of the instrument main plates; 
this is only 1/3 of the temperature range tolerance, +10° F. 
From Figure 6 of Reference 2, the difference in the computed 
from the measured amplitude of node 96 is less than 35 per- 


cent; thus the error in the instrument temperature amplitude 
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is probably on the order of 1° F, and is, in the light of the 


above range tolerance, negligible. 


The Fairing Section 


The fairing or conical section of the "J". system was 
modeled in the same manner described above for the barrel 
section. However, the mathematical model does not include 
the first recovery unit of the system and that unit's influence 
upon the fairing temperature is not accounted for. That unit 
has an ablative, reentry shield that joins the conical skin 
section. The ablative material, a plastic, is not coated with 
the thermal paint mosaic and has a different radiant energy 
absorptivity, 2, and emissivity, ¢€, than that for the rest 
of the vehicle. The value of the ./fe given for a representa- 
tive ablative material in Reference 10, p. F-57, is about 
0.51/0.75 or about 0.68. This value, being less than 1.0, 
will result, according to thé equations of Appendix E, in a 
recovery unit temperature that is below that of the fairing. 
From Vidya Division's thermal library, the radiation heat 
fluxes to a 15° cone for a . = 40° orbit result in the mean 
Orbital hearing fluxes, solar plus albedo, and earthshine of 


G4 G. = 83.7 Btu/hr-ft* ‘ 
G_ = 28.0 ~ 
e 


and, uSing the equation of Appendix E, 


Ce 4 
op = 472° R or 12° F 


“Computations performed by Vidya Division using the 130-node 
model indicate that a 15° F variation in the mean temperature 


of the fairing secticn resuits in an instrument maj 
temperature variaticn @@7Yess than 2° PF. 
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This would represent the orbital mean temperature of the 
external surface of the ablative material. Thus, this ablative 
section can be expected to lower the temperature cf the fairing 
by thermal conduction through the structural elements uniting 
the two sections. The measured temperature for the interior 

of the first recovery unit is about 42° F (Appendix A); the 
thrust cone of this unit has a measured temperature of about 
45° F (Appendix A). These higher interior temperatures result 
from the active heating used in this section and from the heat 
transfer from the fairing section to the first recovery unit. 
Therefore, the measured mean temperatures of the fairing 
section of the vehicle can be expected to be lower than those 
computed. However, the measured transient temperatures are 
functions, also, of the nodal masses as described above for 

the barrel section. The model, having thermal masses that are 
larger than the actuai reactive thermal masses of the system 
may for some assumed thermal environments result in temperature 
predictions that are lower than those measured. Without 
modeling the first recovery unit and the fairing skin structural 
components the algebraic sign of the results compared to that 
for the measurements for the fairing skin cannot be predicted 
since the measured mean temperature, low because of the thermal 
influence of the first recovery unit, may or may not be suf- 
ficiently low compared with the computed value to maintain the 
fairing transient temperature amplitude extremes below the 
maximum or above the minimum values of the computed transient 


4 
temperatures. 


